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Abstract 

O . 

We prove that the r-function of the integrable discrete sine-Gordon model apart 
from the "standard" bilinar identities obeys a number of "non-standard" ones. 
They can be combined into a bivector 3-dimensional difference equation which is 
shown to contain Hirota's difference analogue of the sine-Gordon equation and both 
auxiliary linear problems for it. We observe that this equation is most naturally 
written in terms of the quantum i?-matrix for the XX Z spin chain and looks then 
<^ \ like a relation of the "vertex- face correspondence" type. 

1 Introduction 



Integrable analogues of soliton equations in discrete space-time attract much attention 
for they accumulate developments of the theory of integrable models in the most pure 
form and therefore serve as a good source of new ideas. The activity in the field of par- 
tial difference integrable equations was initiated by Hirota who proposed discrete 
analogues of the KdV, the Toda chain, the sine-Gordon (SG) equations, etc (actually of 
almost all interesting soliton equations). The systematic studies were carried on in the 
works 0, |J, 0, [§]• Among other things, they allowed us to understand why Hirota's 
"empirical" method did really work and generated true discrete analogues of soliton 
equations. 

An important contribution to this understanding is Miwa's interpretation || of the 
discrete time flows, in which discretized integrable equations are treated as members of 
the same infinite hierarchy as their continuous counterparts. The general method to 
produce discrete soliton equations from continuous ones || relies on the Miwa transfor- 
mation - a special change of (in general infinitely many) independent variables connecting 
continuous and discrete time flows. A fundamental new ingredient making all this really 
meaningful is the r-function of an infinite hierarchy of integrable equations ||. Briefly 
speaking, the r-function solves all equations of the hierarchy simultaneously. 
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In few words, the main message of the present paper sounds as follows: without any 
quantization the r-function is already on friendly terms with the quantum i?-matrix, a 
famous object specific for algebraic structure of quantum integrable models. What we 
actually have in mind can be most easily demonstrated for the discrete SG model on 
which we concentrate in the sequel. 

In its original form, Hirota's discrete analogue of the SG equation || reads 

sin (4>(u, v) + 4>{u + 1, v + 1) — 4>{u + — (f)(u, v + 1)) 
= Ksm((f)(u,v) + (j)(u + l,v + l) + (f){u + l,v) + (j){u,v + l)), (1.1) 

where u,v are lattice light cone variables and k is a constant. Recently this equation 



was shown [10 to play a remarkable role in the discrete analogue of continuous two- 



dimensional differential geometry. The Faddeev-Volkov version |]TTJ of this equation, 

ip(u, v + l)ip{u + 1, v + 1) — ip(u, v)ip{u + 1, V ) 
= Kip(u + 1, v)if)(u+ 1, v + 1) — Kip(u, v)if}(u, V + 1), (1.2) 

is sometimes more convenient to deal with. It can be obtained from eq. ( II . 1| ) via the 



substitution ip(u, v) = exp(2ie mv (f)(u, v)). It was shown [[11]], [12| that different continuum 
limits of eq. ( |1.2| ) yield SG and KdV equations. Moreover, eq. ([T2| itself is equivalent to 
Hirota's discrete analogue of the KdV equation (see eq. (|3.9|) below). It should be also 



noted that classical and quantum integrable versions of the SG equation on the space 



lattice with continuous time were suggested in [13| . 

The structure of discretized soliton equations becomes much more transparent if to 
pass to the r-function. From more practical point of view, the reformulation in terms of 
the r-function makes somewhat misteriously looking substitutions easy to guess. 

To be more specific, let us introduce the r-function for eq. ( |1.2p . The idea of Hirota's 
approach is to represent the equation as a consequence of a 3-term bilinear identity for 
the r-function T n (u,v) @], which is simpler and in a sense more fundamental: 

T n (u+l,v)T n+1 (u,v+l) - KT n (u,v+l)r n+ i(u+l,v) = (l-K,)T n (u+l,v+l)r n+1 (u,v). (1.3) 

Here n is a third discrete variable. This equation is valid in the general 2D Toda lattice 
hierarchy |TJ]] discretized along the lines of 0. The subhierarchy with the additional 
constraint 

T n+2 (u,v) = T n (u,v) (1.4) 



is the (discrete) SG hierarchy or, equivalently, 2-periodic 2D Toda lattice [frjj]. Taking 
into account this constraint, it is convenient to write 



r (ii, v) = t(u, v), 7i(u, v) = f(u, v) . (1.5) 

Setting 

// \ r(u,v) 
r(u, v) 

we immediately see that eq. (|1.2|) follows from eq. ( |1.3| ). Higher equations of the hierarchy 



include more discrete variables. 
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In this paper we prove a number of additional identities obeyed by the r-function 
of the SG hierarchy. Among them there are discrete KdV and Toda chain equations 
in bilinear form. Combining them into one matrix identity, we recover that it can be 
most naturally written in terms of the quantum .R-matrix of the XX Z spin chain. The 
identity then looks similarly to what is called "the vertex-face correspondence" between 
models of statistical mechanics on the lattice (see |T5|-P!7|). "Vertex" and "face" types of 



Boltzmann weights are related through " intertwining vectors" . In our case components 
of the "intertwining vectors" are identified with the r-function evaluated at different 
space-time points. 

Treated as an equation, not just an identity, the above mentioned relation contains 
two linear problems and, therefore, implies their compatibility. In its turn, this implies 
a zero curvature condition on the lattice for certain matrix functions L, M (the L-M 
pair), which is equivalent to the discretized SG equation. In this way we reproduce the 



results of our previous papers |]18| , where the L-M pair has been written in terms of the 
quantum i?-matrix. 



2 General form of 3-term bilinear equations 

Here we briefly review some elements of the bilinear formalism for difference soliton 
equations. For more details see e.g. f|, RT9 . 

The key object of the bilinear formalism is r-function. We start with the r-function 
of the 2D Toda lattice hierarchy |p^| . The r-function depends on two infinite sets of 
continuous times t = {t 1 ,t 2 ,t 3 , . . .}, t = {ii, t 2 , t 3 , . . .} (which we call "positive" and 
"negative" respectively) and on a discrete "time" nfZ. For our purpose we do not need 
the "negative" times, so let us write the r-function as r n (t 1; i 2 , t 3 , . . .) = r n (t), where all 
"negative" times (not shown explicitly) are supposed to be fixed. 

First of all we recall how to introduce discrete time flows via Miwa's transformation. 
We use the standard short hand notation 

r(t ± k[z\) = r(ti ± kz, t 2 ± ^kz 2 , t 3 ± ^kz 3 , . . .) , keZ. 

The discrete time k is introduced by the definition 



r, 



n 



r„(t (0) - A^A" 1 ]) , (2.1) 



where stands for constant "background" values of the continuous times and A& G C 
is a parameter of the flow sometimes called Miwa's variable. The A& should be thought 
of as a continuous "label" associated with the discrete flow. From physical point of view 
A^ 1 serves as a kind of lattice spacing in the fc-direction, as it is seen from ( |2.1| ). One 
may introduce arbitrary number of independent discrete flows a,b,c, . . . according to the 
same prescription: 



T a b c... = Tn(t (0) _ a[A -l] _ 6[A -1] _ c[A -l] (2 _ 2) 



It is clear from (|2.2|) that if A a = A&, then the flows a and b must be identified. This 
obvious remark will be important in what follows. 



3 



From now on the discrete times will be written as arguments of the r-function (like 
T n (a, 6, c, . . .)) rather than upper indices. Since the continuous times will not enter 
the play any more, this should not lead to a confusion. For any two discrete flows a, 6 we 
put 

As a function of the discrete times, the r-function obeys a number of bilinear partial 
difference equations. Let us recall their general form. For any triplet {a, b, c} of discrete 
flows one has the following 3-term bilinear equation [Q , |J for the r: 



\b c r(a + 1, b, c)r(a, b + 1, c + 1) 
+ A ca r(a, b + 1, c)r(a + 1, b, c + 1) 

+ X ab r(a,b,c+ l)r(a + 1,6 + l,c) = 0. (2.3) 
For any quadruplet {a, 6, c, d} one has another 3-term bilinear equation: 



A a dAb c r(a + 1, b, c, d + l)r(a, b + 1, c + 1, d) 
+ A M A ca r(a, 6 + 1, c, d + l)r(a + 1, 6, c + 1, d) 

+ A cd A a6 r(a, 6, c + l,d+ l)r(a + 1,6+ l,c,d) = 0. (2.4) 



There are links between eqs. (|2.3j ), ( |2.4|) in both directions. On the one hand, eq. (|2.3| ) 



is a particular case of (O) when A^ — > oo. (According to ( |2.1| ), this limit means that 



the dependence on d in the r disappears.) On the other hand, eq. ( |2.4j ), though linearly 
independent of eqs. (|2.3|), is an algebraic consequense of equations of the type Q2.3p written 
for the triplets {a, 6, c}, {a, 6, <i} and {a, c, <i}. In this sense all what we are going to derive 
in the sequel is implicitly contained in equations of the type ( |2~3| ) . 

Although the discrete time n was distinguished from the very beginning, we can treat 
it on equal footing with the other ones, assigning to it the "label" A n = 0. In particular, 
one can write equations similar to (|2.3|), ( |2.4j ) for the triplet (n, a, 6) and the quadruplet 
(n, a, 6, c) with A n = 0. 



Example. For illustrative purpose we recall the explicit form of the iV-soliton r- 
function. Fix complex parameters Pi,P2, ■ ■ -Pn, Qi, ?2> • • • , Qn an d Q\, ct 2 , • • • , cejy. Then 
the r-function (as a function of discrete variables ai,a 2 ,a 3 , . . .) is given by the determi- 
nant formula lEDfl 



T n (ai, a 2 , ■■■) = det 



Sij + ' 



—) n ~ 



Pj-Qi \<lj) j i\Qj-K 
It obeys all the bilinear relations listed above. 



J l<i,j<N 
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3 Bilinear equations for the r-function of the dis- 
crete SG model 



Let us specify the general formulas for the case of the discrete SG model. Fix two 
discrete flows u, v with the corresponding parameters X u , X v . For brevity we set // = X u , 
v = X v . Let T n (u, v) be the r-function. As it was already said in the Introduction, the 
SG hierarchy is defined by the reduction condition ( |1.4|) . (In the example given at the 
end of the previous section this reduction is equivalent to the constraints Pi + qi = 0, 



Using the notation ( |1.5| ), we rewrite eq. (|2.3|) for the triplet (n, u, v) as the pair of 
coupled equations for r and f: 

ut{u + 1, v)t(u, v + 1) — /if (u + 1, v)t(u, v + 1) = (y — n)t{u, v)r(u + 1, v + 1) , (3.1) 

uf{u + 1, v)t(u, v + 1) — jiriu + 1, v )f(u, v + 1) = {y — /u)r(w, v)f{u + 1, v + 1) , (3.2) 
They can be rearranged to the following equivalent form: 

fir (u, v)t(u + 1, v + 1) + z/t(m, v)f{u + 1, t> + 1) = [y + h)t(u, v + l)f(u + l,v) , (3.3) 

yur(u, t>)f(u + 1, v + 1) + z/f(u, v)t{u + 1, f + 1) = [y + /i)f(u, v + l)r{u + l,v) . (3.4) 

The Faddeev-Volkov equation ( |1.2|) for ip{u,v) = t(u,v)/t(u,v) holds with k = fi/u. As 
is already clear from the previous section, the parameters \i and v play the role of inverse 
lattice spacings in the chiral directions u, v respectively. 

There is a relation connecting values of r and f along the w-direction (or the in- 
direction). To derive it, introduce a third discrete variable u' with Miwa's variable //. 
Then equations of the form (|3.3|), ( |3.4| ) are valid for u' in place of v (and // in place of 
fi) while v is not shifted and enters as a parameter. Recalling the remark after eq. ( |2.2| ) . 
we should identify u' with u as // = /i; this yields the relation 



We have obtained basic bilinear equations of the discrete SG model. 

Besides the basic equations given above, the r-function obeys some other 3-term 
bilinear identities which follow from the basic ones. Here we mention three of them: 



i = 



1,...,N.) 



t{u — 1, v)t(u + 1, v) + r{u — \,v)t{u + 1, v) — 2t(u, v)f{u, v) . 



(3.5) 



(/i + z/)r(w — 1, v + 1)t{u + 1, v ) + (/i — u)t{u + 1, v + 1)t(u — 1, v ) 
2ht(u, v)t(u, v + 1) , 



(3.6) 



(/i + z/)r(u — 1, v + l)f{u + 1, v) — (fj, — u)t{u + 1, v + l)f(u — 1, v) 
2ur(u, v)f(u, v + 1) , 



(3.7) 
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(fi+u) 2 T(u-l,v + l)r(u + l,v-l) - (h-v) 2 t(u+1,v+1)t(u-1,v-1) 
= ^u(t(u,v)) 2 . (3.8) 



Their derivation from the basic equations is given in the Appendix. Note also that since 
r and f differ only by the shift of the hidden variable n and r„ +2 = r n , the exchange 
r <-> f in eqs. fl3.6|) - ([3.8|) automatically yields true identities which we do not write down 
explicitly. 

We note in passing that eqs. ( |3.6|) and (|3.8|) provide bilinearization of Hirota's discrete 
KdV equation and ID Toda chain in discrete time respectively. Indeed, it follows from 
that 

r(«,»M«+l,«, + l) 

T{U + 1, V)T[U, V + 1) 

solves the discrete KdV equation |L| 

U(u+ l,v + 1) - U(u,v) = l^(u-\u,v + 1) - U-\u + l,v)) (3.9) 

1 — K V ' 

with k = fi/u. Note also that substituting 

+ 1) - «V>(u + l,v) 

U (U, V) = — r , 

(1 — K)1p{U, V) 

we get the Faddeev-Volkov equation ( |1.2|) for ip. 
Further, eq. (|3.8|) implies that 

r(u + l,u + l) 
¥?(u,i;) = log r 

T{U, V) 

obeys the equation 

exp(tp(u + 1, v — 1) + (^(^ — 1, f + 1) — 2</?(w, w)^) 
1 + g exp((p(u + 1, v + 1) — u)) 



1 + (?exp((^(-u, t>) — </)(m — 1, v — 1)) 



(3.10) 



with g = Afiu/i^fi — z/) 2 , which after passing to the "laboratory" coordinates |(u ± v) 
coincides with the ID Toda chain equation in discrete time |2| (see also |yj). We stress 
that all this holds for the same r-function t(u, v) of the discrete SG equation. 



4 3D difference equations for the SG r-function 

Let us turn to bilinear equations for the SG r-function which include more discrete 
variables. In some sense they are discrete counterparts of higher equations of the SG 
hierarchy. Here we consider the case of one such "extra variable" w introduced according 
to (|2.2| ). Set \ w = (. The r-function is now written as t (u,v,w) = r(u,v,w) and 
Ti(u,v,w) = f(u,v,w). It satisfies the equations 
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(C - h)t(u, v, w)t(u + 1,v,w + 1) 
= (t(u + 1, v, w)t(u, v, w + 1) (4.1) 

— [ir{u + 1, v, w)t(u, V, w + 1) , 

(C — ^)r(u, u, w)t(u, v + l,w + 1) 
= (t(u,v + l,w)f(u,v,w + 1) (4.2) 

— z/f(w, t> + 1, w)t(u, V, w + 1) 

(and the same with the change r «-> f). They are nothing else than eq. (|3.1| ) written for 



the triplets (p,u,w) } (p,v,w) respectively. Similarly, the analogues of eqs. fl3.3|) , ( |3.4| ) are 



(C + yu)r(-u + 1, v, w)f(u, v, w + 1) 
= /if(ii + 1, t>, w + l)r(u, f , w) (4.3) 
+ (t(u + 1,v,w+ l)f(u,v,w) , 

(C + v)t(u, v + 1, w)f(u, f , w + 1) 
= z/f(w, t> + 1, w + 1)t(w, t>, u?) (4.4) 
+ (t(u, v + 1, w + l)f(u, v, w) 



(and again the same with the change r <-> f) 
Eq. (|2.3j ) for the triplet (w, t>,u>) gives 



(C - p)r(u, v + 1, u;)r(w + 1, t> , u; + 1) 
— (C — i/)r(w + 1, v , w)t(u, v + 1, w + 1) (4.5) 
+ (/i — ^)r(u, u, w + l)r(u + 1, t> + 1, u>) = . 



Eq. (|2.4p for the quadruplet (p,u,v,w) gives 



z/(C — fi)f(u, v + 1, w)t{u + 1, v, w + 1) 
— /^(C — ^)f(w + 1, i>, w)r(u, v + 1, u> + 1) (4.6) 
+ ~~ u )t(u, v,w + 1)t{u + 1, v + 1, u>) = . 

In the next section we need two more complicated bilinear equations: 

(/i + v){C - fi)r(u — l,v + l, w)t(u + 1, v, w + 1) 

+ ( / u-z/)(C + / u)r(w- 1 ^ ) ^ + 1 ) r (w + 1 ^ + 1 ) « ; ) ( 4 -7) 
= 2//(C — v)t{u, v, w)t(u, v + 1, it; + 1) , 



(fi + v)(C - p)t(u — 1,v + 1,w)t(u + l,v,w + 1) 
(;U - u)(( + /i)r(u - w + 1)t(u + l,v + l,w) 

t(m, f , w)T{u } v + l,w + 1) 



r(u, v, w)f(u, v + 1, w + 1) 
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(and the same with r <-> f). They can be derived from eqs. ( |4.1|) -(0) in the manner 
similar to how eqs. fl3.6|) - fl3T8"D have been derived from the basic equations (|3.1|) - (|3^) . 
The derivation is outlined in the Appendix. 



5 i?-matrix form of the 3D bilinear identities 



Our goal in this section is to combine the identities (|4.7| ), ( f4.8| ) into one compact equation 
suggesting parallels with lattice integrable models of statistical mechanics. 



Let us introduce the matrix of Boltzmann weights of the 6- vertex model: 



R 



( a ^ 

b c 

c b 

V a / 



and call it i?-matrix. At the moment we do not imply any specific parametrization of 
the matrix elements. Assuming the standard block structure of the i?-matrix, we say 
that it acts in the tensor product V\ ® V 2 of two 2-dimensional spaces isomorphic to C 2 : 
R G End(y x g> V 2 ), V 1 = V 2 = C 2 . In other words, 

R = ^(a + b)J ® I + i(a - b)a 3 <g> a 3 + ^ c (°"i ® <J\ + o 2 ® o 2 ) ■ 
(Here I is the unit matrix and Oi are Pauli matrices.) Hence the i?-matrix acts to tensor 



product of the vectors X 



X 2 



RX®Y = R 



eV u Y 

( xiVi \ 
x 2 yi 
X1V2 
V x 2 y 2 ) 



yi 

V2 



G Vo as follows: 



( arcij/i \ 
bx 2 yi + cxiy 2 
cx 2 yi + bxiy 2 

V a^22/2 / 



It is easy to verify that eqs. (|4.7|), ( f4.8[ ) are equivalent to the relation 
t(u,v,w) \ I t(u, v + 1, w + 1) 

Rl ® 

r(u, 1>,iu) / y — f(tt, f + 1, w + 1) 

■(« — 1, v + 1, iu) \ / r(w + 1, v, w + 1) 
<\ I ® 

f(ii — 1, f + 1, w) / \ — f(lt + 1, V, w + 1) 



+ 



I t{u + 1, v + 1, w) 
\ —t(u + 1, v + 1, w) 



I r{u — 1, v, w + 1) 
\ t(m — 1, w, w; + 1) 



(5.2) 
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where the i?-matrix has elements 



a = 2MC 2 - ^) 
b = 2z/(C 2 - fi 2 ) 
c = 2C(^ 2 - z/ 2 ) 



(5.3) 



and a, [3 are scalar coefficients: 



I(a + b + c) = ( M + ^)^ + ()((- )U ) 



/? = i (a - b - c) = - i/) (1/ + C) (C + A*) 



(5.4) 



After a proper rescaling and extracting a common multiplier, R becomes the standard 
trigonometric i?-matrix with spectral parameter £: R = 2fiu(R K ((/fi), where k = fi/v 
and 

/ k^-kT 1 ^ 1 \ 

2 — Z^ 1 K — K 

K-K- 1 Z-Z' 1 

kz-k~ 1 z~ 1 ) 

is the i?-matrix of the XXZ-type with the quantum deformation parameter k and the 
spectral parameter z. It obeys the Yang-Baxter equation 



R K (z) 



V 



(5.5) 



(5.6) 



where R^ acts as R K in the tensor product of z-th and j-th spaces and as the unit matrix 
in the remaining space. 

To clarify the meaning of eq. ( |5.2|) , we rewrite it as follows: 

R K (z)$(u,v) <g> ^(u,v + 1) 
= a(z)$(u-l,v+l) <g) ty(u+l,v) +p(z)<T3$(u+l,v+l) ®a 3 ty(u-l,v) . (5.7) 



Here 



/ t(u, v,w) 



$(u,v) 



(w is supposed to be constant), z = and 

a(z) = |(1 + k- 1 *- 1 )^ + l)(z-l), 

p{ z ) = Ki + ft- 1 ^- 1 )^- 1)(^ + 1). 



r(w, v, w + 1) 
-f(u, u, iu + 1) 



(5.8) 



(5.9) 



Note that \l/ implicitly depends on z through the shift of w. In this form the identity 
resembles relations of the "vertex-face correspondence" type. It contains two linear prob- 
lems whose compatibility condition is the discrete SG equation with \1/ playing the role 
of the " wave function" . 
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ty(u,v + 1) 




u 



ty(u + l.v) 



To see this, let us take the scalar product of both sides of eq. (|577| ) with the vector 
cti$(m + 1, f + 1) in the first space, orthogonal to a^{u + 1, v + 1). We get: 



(^(u + M + l), R k (z)$(m,u)) i 
oi$(u + l,f + 1), $(« - l,v + 1) 



*(«, w + 1) = + 1, f ) 



(5.10) 



The subscript 1 in the numerator means that the scalar product is taken in the first space 
V\, so the result is a 2x2 matrix in the second space. In a similar way, taking the scalar 
product with the vector — ia 2 ^{u — l,v + 1), we get: 



(<7 2 $(u-l,i> + l), R K (z)$(u,v)) i 

+ 1, V + 1), - 1, V + 1) 

Note that 



*(u, f + 1) = z/3(z)<7 3 *(tt - 1, u) . 



(5-11) 



o"i$(u + l,f), $(m - l,ti)) = 2t(u,v 



)r(u,v) 



due to ( |3.5|) . Eqs. ( |5. 10| ) , (|5.11|) are two linear problems for the transport of the wave 
function \1/ in the two directions displayed in the figure. The left hand sides of these 
equations provide us with the L-M-pair depending on the spectral parameter z that 
comes from the quantum i?-matrix. At the same time this spectral parameter coincides 
(up to a constant scale factor) with Miwa's variable (. The compatibility of the two linear 
problems is expressed as a discrete zero curvature condition equivalent to the discrete 
SG equation. Given a solution to (|5.2|), at any slice w = const we have a solution to 
the discrete SG equation in bilinear form. Therefore, we have got another proof of the 
results of [IS], where local L and M matrices for discrete time shifts were expressed in 
terms of the quantum i?-matrixQ. 



1 Comparing with the final formulas of |iq |, there is a minor difference since here we actually deal 
with L^ 1 rather than L. 
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6 Discussion 



We have shown that the r-function of the discrete SG hierarchy apart from the standard 
bilinear equation (|1.3f) satisfies additional bilinear relations ( ^6|) -( ^78|) , the first (respec- 



tively, the last) of which provides bilinearization for the discrete KdV equation (respec- 
tively, the ID Toda chain in discrete time). Besides, we have proved some additional 
bilinear relations for the "extended" SG r-function depending on 3 discrete variables. 
They are combined into the matrix relation (|5.7|) which contains the quantum i?-matrix 
of XXZ-type and resembles the vertex-face correspondence. In less dry terms, this fact 
would probably mean that the structure of integrable bilinear partial difference equations 
is to high extent determined by solutions of the quantum Yang-Baxter (YB) equation 
T6I). 



Recall that both classical and quantum YB equations already appeared in connection 
with integrable mappings. First, the (classical) r-matrices are known to play a signifi- 
cant role in the algebraic structure of discrete time integrable mappings |21| , |[22|| . The 



quantum YB equation as arising in purely classical problems was discussed in the works 
f23fl , pi]] , though the class of relevant solutions seemed to be very far from i?-matrices of 
known quantum integrable models. In the present paper we have observed that a well 
known particular solution to the quantum YB equation is really relevant to the discrete 
SG hierarchy although in a different sense. However, in our construction the role of the 
quantum YB equation itself is not yet clear enough. Hypothetically, it is responsible for 
the fact that the linear problems ( |5.10 ) and ( 5.1 1| ) are not too overdetermined. (In other 



words, the equations obtained as their consistency conditions do not contradict to each 
other.) 

At last, we should stress that the "quantum parameter" k of the trigonometric R- 
matrix in our context seems to have nothing to do with any kind of quantization. It is 
related to the mass parameter and the lattice spacing of the classical model. 
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Appendix 



Here we show how to derive the identities ( |3.6] )-(|3~8|) and ( [4.7| ), 



Consider eqs. (|3.1|) , ( |3.3j ) with the shift u — > u — 1 in the latter: 

vr{u + 1, v)f(u, v + 1) — fif(u + 1, v)t(u, v + 1) = {y — fi)r(u + 1, v + l)f(u, v) 

vr{u — 1, v)f(u, v + 1) + fif(u — 1,v)t(u, v + 1) = {v + fi)r(u — 1, v + l)f(u, v) 
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Rewriting them system of linear equations, 
vt(u+1,v) —fj,f(u+l,v) \ / t(u,v+1) 

vt(u—1,v) /j,t(u — 1,v) J \r(u,v+l) 
and using the Kramer rule, we get: 



--f(u, v) 



(i/—fi)r(u+l,v+l) 
{v+ijl)t{u—1,v+1) 



f(u, v + 1)D = f(u, v) 



t(u, v + 1)D = f(u, v) 



(u— //)r(w+l, v+1) —/jlt(u+1,v) 

(v+/jl)t(u—1,v+1) /if(u—l,v) 
vt(u + 1,v) (u — fj)r(u+l,v+l) 
ut(u—1,v) (v+h)t(u — 1, v + 1) 



(Al) 



(A2) 



where 



D = v/jl(t(u+1 : v)f(u — l,v) + T(u—l,v)f(u+l,v) \ . 

Taking into account eq. (|3.5|) , we see that D = 2ufj,r(u, v)f(u, v) and, therefore, and 
(|A2|) coincide with (|3.7|) and ( |3 .6|) respectively. Eq. ( p.8|) is obtained in a similar way 
from eq. ( p.7| ) and its counterpart with exchanged flows u and v. 

To prove ( fOD , consider eq. (|3.7|) and its counterpart with v, v substituted by w, (: 



(C + h)t{u — 1, v, w + 1)t(u + 1, v, w) + (C — h)t{u + 1, v, w + 1)t(u — 1, v, w) 
= 2(t(u, v, w)t(u, v, w + 1) . (A3) 
These two can be rewritten as the following system of linear equations: 



/ (u+ijl)t(u— 1, v+1, w) {y— /i)t(u+1, v+1, w) 



I t(u+1,v,w) 



\ ((+ii)t(u— 1, v, w+1) ((—ii)t(u+1,v,w+1) J \ t(u—1,v,w) 
I vt(u,v+1,w) 

2t(u, v, w) 

V (r(u,v,w+l) 



(A4) 



Solving this system for f(u + 1, v , w) and taking into account eq. (4.6), we get 

To prove ( |4.8| ) , consider eq. ( |3.(j| ) with r — > f and eq. ( |A3[ ) combined into the following 

system: 

/ {v+h)t(u—1,v+1,w) (//— v)t(u+1, v+1, w) \ I t(u+1,v,w) 

V (C+A t ) r ( M -l)' u ) w + 1) ((—h)t(u+1,v,w+1) J \ t(u—1,v,w) 

I /j,f(u,v,w)f(u,v+l,w) \ 
= 2 . (A5) 

\ (,t(u, v, w)f(u, V, w + 1) J 

Proceeding similarly to the previous case and using equations from Sect. 4, we come to 
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